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We show how the continuity equations expressing conser- 
vation of topological point or string defect charge can be used 
to determine the order-parameter correlation function for the 
phase-ordering kinetics of the 0{n) model in the special case 
where the order parameter is constrained to be near a defect 
core. In this regime we find a self-consistent solution by as- 
suming the order parameter is Gaussian. The resulting linear 
equation for the order-parameter correlation function has as 
its solution the Ohta-Jasnow-Kawasaki form. 

PACS numbers: 05.70.Ln, 64.60. Cn, 64.75. -f-g, 98.80. Cq 



In this letter we show that the phase-ordering kinetics 
of defects in the 0{n) model can be self-consistently eval- 
uated assuming the order parameter is a Gaussian field 
near the defect cores. This, in turn, clarifies the nature of 
recent determinations of defect correlation functions for 
the 0{n) model We consider systems with either 

point defects (n = d, where d is the spatial dimensional- 
ity) or string defects {n ~ d—l). In addition to their im- 
portance in condensed matter, these systems are also rel- 
evant to problems in cosmological structure formation . 
Our approach is based on the recently derived continu- 
ity equation expressing conservation of topological defect 
charge and differs from previous calculations which make 
more direct use of the time-dependent Ginzburg-Landau 
(TDGL) equation for the non-conserved order parame- 
ter [^,0 . The continuity equation for systems supporting 
string defects has not appeared previously and is estab- 
lished here. These continuity equations are used to find 
the equation satisfied by the order-parameter correlation 
function where the order parameter is restricted to be 
near a defect core. Since it is easy to confuse the usual 
order-parameter correlation function with the correlation 
function for the order parameter constrained to be near 
a defect core, we will, for reasons which will become clear 
below, refer to this constrained quantity as the auxiliary 
field correlation function /. We obtain a linear equation 
for / which has as its solution the Ohta-Jasnow-Kawasaki 
(OJK) 1^ form. Since the equation is linear, the Gaus- 
sian assumption is consistent. 

Most of the focus in phase-ordering kinetics has been 
on developing theories for the order-parameter correla- 
tion function. The belief has developed that we have a 
fairly good understanding of how to calculate this quan- 
tity in the late-time scaling regime. Two theories have 
evolved which map the order parameter field onto an aux- 
iliary field and invoke the so-called Gaussian closure ap- 



proximation. We shall refer to the correlation function 
for this auxiliary field as fop- The OJK method, with 
elaborations by others leads to a simple form for 

the auxiliary field correlation function. The second ap- 
proach, originated by Mazenko determines fop self- 
consistently through the solution of a nonlinear eigen- 
value problem. A key point we make here is that it is 
not necessary that / = fop- 

Since the order-parameter correlation function is a 
rather structureless quantity which does not give a great 
deal of direct information about the underlying disor- 
dering agents, Liu and Mazenko Q examined the cor- 
relations between the defects themselves. The key ele- 
ment in this work was, following Halperin [ pl| , to iden- 
tify the positions of the defects with the zeros of the 
order-parameter field which could, in turn, be mapped 
onto the zeros of a Gaussian auxiliary field. The point 
defect charge density correlation function Gp was deter- 
mined using this technique in terms of /, assuming that 
/ — fop- For large defect separations, Gp was found to 
agree with numerical simulations p2| and experiments 
||l3| relevant to the case n = 2. However, for n = 2 
the theory produced an unphysical diverg ence in Gp at 
short-scaled distances x due to a non-analytic piece in 
auxiliary field correlation function fop at small x. This 
non-analytic piece occurs in the most elementary self- 
consistent theory for /op but Mazenko and Wickham 
§ have since shown that one can construct the theory so 
that fop is analytic in x. This development highlights 
the difference between the order-parameter correlation 
function, whose short distance non-analyticities are es- 
sential and lead to the observed generalized Porod's law 
for the structure factor, and / which must be smooth 
in order to have sensible theories of defect position and 
velocity correlations. The / obtained from OJK is 
smooth. 

The approach we take in this paper is distinct and in- 
dependent from that used by OJK or that developed in 
1^. The technique used here is directly based on con- 
siderations of topological charge conservation and is less 
directly based on the TDGL model. This leads to the im- 
portant difference that we treat quantities that are con- 
strained to be evaluated at the defect core. Although 
we recover the OJK result we obtain information on the 
regime of applicability of this result. 

The system studied here has a defect dynamics gen- 
erated by the TDGL model for a non-conserved n- 
component vector order parameter 'ip{f,t): 
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dt Sip 



(1) 



where F is a kinetic coefficient, is a Ginzburg-Landau 
effective free energy assumed to be of tlie form 



F 



(2) 



Ttie coefficient c is positive and the potential V is as- 
sumed to be of the 0(n)-symmetric, degenerate double- 
weU form. We assume that the quench is from an initial 
high-temperature disordered state to zero temperature 
so the usual noise term in (|l]) is set to zero. 

In previous work j^] on the order-parameter correla- 
tion function progress was made by mapping the order 
parameter ip onto an auxiliary field m, with the require- 
ment that away from the defect cores 



(3) 



where ipQ is the magnitude of ■0 in the ordered phase. 
Physically, we interpret m to be the position relative to 
the nearest defect and expect that near the defect cores 



ip — am + brfi[Tnf' + 



(4) 



where the coefficients a and h depend on the details of 
the potential V. Equations (^) and (P represent topo- 
logical charge ±1 defects, which have the lowest energy 
and dominate the late-time regime p^ . In the theory 
for order-parameter correlations j7| property (|^) is cru- 
cial, whereas, in the theory of defect motion presented 
here property (^) is relevant since we always work near 
the defect cores. 

In the simplest models f^J^, considered here, the aux- 
iliary field rfi is assumed to be a Gaussian field with a 
normalized correlation function / defined as 



(m^(l)m„(2)) 
V^o(l)^o(2) 



(5) 



with (5p,yS'o(l) — (mp(l)my(l)). Here we use the short- 
hand ?T7,^((1) — m^((ri,ti). It is well-established that for 
late times t following the quench the dynamics obey scal- 
ing and the system can be described in terms of a single 
growing length which is characteristic of the spacing 
between defects. For a non-conserved order-parameter 
Ht) ~ t^/^ at late times. Since we interpret to be 
the distance to the nearest defect we expect \rh\ ~ L. 
At equal times (ii = t2 = t) in the scaling regime the 
auxiliary field correlation function can be written solely 
in terms of the scaled length a; = |r2 — ri|/L(i). Hence 

The emphasis in this letter is on defect densities like 
the charge density for point defects, given in terms of the 
order parameter by 



(6) 



where the Jacobian associated with the change of vari- 
ables from the set of defect positions to the field %p is 
defined by 

6/^1/^2 • Mil i^ the n-dimensional fully anti-symmetric tensor 
and summation over repeated indices in (|^) is implied. 

It was shown in a direct manner in Q that p satisfies 
the continuity equation for topological charge 



where the current is defined as 



(8) 



J, 



1 



(9) 



The derivation of (^) is independent of the details of the 
TDGL model, except that equation ([^) is first order in 
time. Since j'^^ is multiplied by the defect-locating 5- 
function we can replace K in J^^^ by the part of K which 
does not vanish as '0 ^ 0. For a non-conserved order 
parameter this means that we can set K = TcV'^ip in 
. Equation (||) is in the standard form of a continuity 
equation, allowing us to identify the vortex velocity field 
as 



V 



(10) 



where it is assumed that the velocity field is used inside 
expressions multiplied by the vortex-locating 5-function. 
For string defects the charge density is a vector given 

by 



Pa = 6(lp)ljJo 



(11) 



with 



(12) 



As in one can obtain the continuity equation satisfied 
by Pa by combining the two identities 



and 



to obtain 



K^UJa = J^^^P'ip'j 



dpa 

dt 



(13) 
(14) 

(15) 
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The string defect current tensor J^^' is defined as 



j: 



al3 



{n-iy. 



(16) 



For the important case of n — 2, d = 3 we can write 

and identify the string velocity as 
1 



(17) 



(18) 



where g = e,.! ^2X^1 Vi/';.2 0- 

We use the continuity equation for topological charge 
to determine the auxiliary field correlation function. We 
begin by examining point defects and require that the 
exact equation 



d_ 

dt 



{p{l)p{2)) ^vf{5m)]Jl->{l)p{2)) 



{K), 



be satisfied at equal times. The presence of the 5- 
functions in (|l^) enables us to use relation (0) to replace 
tj} with the Gaussian auxiliary field to in (|l9|). The left- 
hand side of (^9|) involves the point defect charge density 
correlation function defined by 



G,(12) - (p(l)p(2)). 



(20) 



As shown in Q], Gp factors into a product of Gaussian 
averages which can be evaluated using standard methods. 
In the scaling regime Gp has the form 



Gp(12) 

with g(x) given by 

g{x) = n\ 



I2n 



h{x) 



dh{x) 
dx 



We define 



h = 



27r 



(21) 



(22) 



(23) 



with 7 = 1/ ■\/l — /2 . It is therefore clear that the left- 
hand side of ( p^ ) is a complicated non-linear function of 
/ and its derivatives. The right-hand side of (|9|) involves 
the evaluation of 



Np{l2)^{5[m{l)]V{l)5[m{2)]jf\2)) 



(24) 



which, like the evaluation of Gp, factors into products of 
averages over the n separate components of to. We easily 
find 



iV^(12) = n!rcB(A)"-i% 
where, in the scaling regime, A is given by 

1 h 



and B is 



A = 



B 



L2 



B 

Z3 



with 



2nB = 



dx 
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(25) 



(26) 



(27) 



(28) 



We have defined S^'^'> = J^afiii'^ a'mp]'^) /n^ and a = 
Sq/L'^, which are both constants at late times. 

Having compiled these results, it is easy to see that 
dlOl ) reduces, in the scaling regime, to 



p— f/i" + a;/i"-i/i'] 

L J 



dx ^ dx 

where we have defined the constant 



P : 



LL 
2fc' 



Equation (p9|) can be integrated to give 



p—ixh) — B 
dx 



(29) 



(30) 



(31) 



with the integration constant determined to be zero by 
the condition that / and h vanish as a; ^ cxd. Equation 
( ^ ) shows us that B is the derivative of a quantity which 
vanishes at a; ^ 00. Thus we can integrate (p3l|) yet again, 



and use (23), to obtain the remarkable final result 



- pxf = vV 



-/ 



(32) 



which is linear in /. 

For string defects, the analogous calculation involves 
evaluating the averages in 



d_ 

m 



{Po.[i)Pfm = vW(j[v^(i)] jW(i)p^(2)) 



+v^^Hpo.{mm]jf^\2)). 

The string charge density correlation function 
G„/5(12) = {pa.{l)pp{2)) 



(33) 



(34) 



was first worked out by Liu and Mazenko j^] with the 
scaling result 

Gaf}{l2) [Gt{x) [Sap - XaXp) + GL{x)XaXp] 



(35) 
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where the transverse function is 
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Gt{x) 



dh 
dx 



and the longitudinal function is 
Gl{x) =n\'^ 



(36) 



(37) 



h is defined in (p3|). We evaluate the average appearing 
on the right-hand side of ( ^3| ) using the same techniques 
as used in the point defect case and easily obtain 



(38) 



where A is defined i n (P6|) and B in (27). Substitution of 
the results (|3^) and (|38[) into ( |33[ ) leads to separate equa- 
tions for the longitudinal and transverse components. 
Both equations reduce to (^). Thus the same linear 
equation determines / for both point and string defects. 
The solution of (Is^) is of the OJK form 



(39) 



/(x) = exp--a; 
where we have used the relation 



(40) 



The constant ji is fixed through a choice of the length 
scale L, since L = y/W^t^^^. We are therefore lead di- 
rectly, without any approximation except that the order 
parameter field can be treated as Gaussian near its zeros, 
to a self-consistent result. 

It should be noted that at leading order in a systematic 
large- iV approximation scheme for a colour index TV Bray 
and Humayun pO| ] were able to recover the previously 
ad hoc Oono-Puri ||^ extension of OJK which, in turn, 
gives exactly (|3^). In the development here, equation (^) 
arises without the need for any such approximations. 

This work highlights the fact that the original OJK re- 
sult is implicitly derived near the defect cores It also 
emphasizes the distinction between theories of the order 
parameter correlation function, such as in Ref. ||^, and 
theories concerning defect correlation functions, where it 
appears one can self-consistently use the Gaussian closure 
approximation if one uses the OJK form for the auxiliary 
field correlation function. The theory of order-parameter 
correlations does a superior job determining the non- 
equilibrium exponent governing the decay of two-time 
autocorrelation functions but the use of the the- 

ory presented here for defect correlations, with its smooth 
O JK-like auxiliary field, avoids the difhculties due to non- 
analyticities at small x. 
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